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A linear open quantum system consisting of a harmonic oscillator linearly coupled to an infinite 
set of independent harmonic oscillators is considered; these oscillators have a general spectral density 
function and are initially in a Gaussian state. Using the influence functional formalism a formal 
Langevin equation can be introduced to describe the system's fully quantum properties even beyond 
the semiclassical regime. It is shown that the reduced Wigner function for the system is exactly 
the formal distribution function resulting from averaging both over the initial conditions and the 
stochastic source of the formal Langevin equation. The master equation for the reduced density 
matrix is then obtained in the same way a Fokker-Planck equation can always be derived from a 
Langevin equation characterizing a stochastic process. We also show that a subclass of quantum 
correlation functions for the system can be deduced within the stochastic description provided by 
the Langevin equation. It is emphasized that when the system is not Markovian more information 
can be extracted from the Langevin equation than from the master equation. 

I. INTRODUCTION 

Feynman observed long ago that the dynamics of an open quantum system may be described in terms of an equivalent 
stochastic problem [Q. In this paper, we elaborate on this insight by showing that a certain class of quantum Green 
functions may be obtained directly as ensemble averages in the stochastic formulation. Unlike earlier treatments of 
the origin of stochasticity in the semiclassical limit of quantum theories [|| we do not assume decoherence; rather, we 
are concerned with the full quantum dynamics of the system. 

The stochastic treatment follows from the observation that the Wigner function of the open quantum system may 
be represented as the formal distribution function resulting from averaging the solutions of an appropriate Langevin 
equation both over the initial conditions and the stochastic source. This representation, which is exact for linear 
systems, may be extended to nonlinear cases through pcrturbative methods 

We compare the main approaches to the analysis of the stochastic dynamics, namely the Langevin and Fokker- 
Planck equations, to the corresponding quantum approaches, namely the master equation and the Wigner function. 
We show that the Fokker-Planck equation is the transport equation for the full Wigner function, and as such it is 
equivalent to the master equation. The Langevin equation, on the other hand, provides a more detailed description of 
the dynamics, in the sense that the class of correlation functions which may be retrieved from the Langevin equation 
is larger than the corresponding class for the master or Fokker-Planck equations unless the dynamics is Markovian. 

The relationship between these different approaches is summarized in the diagram of Fig. 1. The most fundamental 
description of the open quantum system is that provided by the Feynman- Vernon influence functional From this, 
we may derive the master equation, which gives the dynamics of the reduced density matrix. The integral transform 
linking the reduced density matrix to the reduced Wigner function allows us to convert the master equation for the 
former into a transport or Fokker-Planck equation for the latter. Feynman's insight that the influence functional may 
be thought as well as an ensemble average over an equivalent stochastic noise allows us to retrieve some correlation 
functions of the quantum problem directly in terms of stochastic averages, while the Fokker-Planck equation associated 
to the corresponding Langevin equation gives back the master equation. Thus, the Langevin equation is a very useful 
tool to gain information on the quantum properties of the system even beyond the semiclassical regime, when it no 
longer describes the actual trajectories of the system. 



*Also at Institut de Fisica d'Altes Energies (IFAE), Barcelona, Spain. 
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FIG. 1. Diagram showing the interconnections between different quantum properties of an open quantum system on the one 
hand, and between the elments of the stochastic description on the other hand, as weU as the connection between both levels 
of description. Quantum objects (denoted by label Q) are placed at the top of the diagram whereas those objects associated 
to the stochastic description (denoted by label S) are placed at the bottom. The equivalence of the reduced Wigner function 
with a distribution function in the framework of the stochastic description provided by the Langevin equation reflects the 
connection formally established between both levels of description. One can take advantage of this connection by working 
within the stochastic framework to obtain relevant information about quantum properties. In the stochastic context, one can 
always obtain a Fokker-PIanck equation for the distribution function associated to the stochastic process resulting from a given 
Langevin equation. From the connection previously stated, the Fokker-PIanck equation can be translated into the master 
equation for the reduced Wigner function or, equivalently, for the reduced density matrix. Similarly, correlation functions for 
the stochastic processes which are solutions of the Langevin equation are closely related to quantum correlation functions for 
system observables. 



The plan of the paper with a summary of the relevant results is the following. In Sec. ^ we briefly review the 
quantum Brownian motion (QBM) model describing a harmonic oscillator coupled bilinearly to a bath of infinite 
harmonic oscillators initially in a Gaussian state and with an arbitrary spectral distribution (general environment). 
We summarize the main formulas that beginning with the infiuence functional lead to the reduced density matrix 
operator and to the master equation. We also give the equivalent evolution equation for the reduced Wigner function. 
These results are well known and the computational details, which can be found in the references provided, are 
omitted. 

In Sec. [II we derive one of the main results of this paper. We show that the reduced Wigner function can be 
written as a formal distribution function for the system variables defined as follows. We consider the ensemble of 
systems in phase space which obey the formal Langevin equation for a given realization of the stochastic force and 
arbitrary initial conditions; the average both over the initial conditions (which involves the initial reduced Wigner 
function) and the different realizations of the stochastic source gives just the reduced Wigner function. The key 
technical point that makes this result relatively easy to derive is the computation of the path integral defining the 
reduced Wigner function in terms of variables which include the initial conditions and the stochastic force; for this 
we need to introduce a functional change which involves the formal Langevin equation. The details of this derivation 
can be found in Appendix As far as we know this result is new. We are aware of a related result by Halliwell and 
Zoupas Q in the limit of large times; see also Ref. [||. 

Having shown that the reduced Wigner function is a formal distribution function, it is clear that the dynamical equa- 
tion for the reduced Wigner function can be deduced using the usual techniques employed to derive the Fokker-PIanck 
equation Q . We recall that whereas a Langevin equation describes how an individual system evolves stochastically, a 
Fokker-PIanck equation describes how the distribution for an ensemble of systems evolves deterministically in phase 
space. The details of this derivation are left to Appendix This constitutes an alternative derivation of the master 
equation to those given by Hu, Paz and Zhang Q (see also ReL 1^), and Halliwell and Yu We also note that this 
process is not reversible in the sense that there may be many Langevin equations that give the same master equation. 
This just reflects the fact that when the dynamics is not Markovian (in the sense of Ref. |]lO|), more information can 
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be extracted from the Langevin equation than from the master equation. 

In Sec. rV we obtain the second main result of this paper. We show that quantum correlation functions for the system 
variables can be obtained within the stochastic description provided by the Langevin equation by explicitly computing 
the closed time path (CTP) generating functional for the system. It turns out that this generating functional can be 
written as an average over the initial conditions times a term that depends on the noise kernel, which contains the 
information on the fluctuations induced by the environment on the system. 

We also show that quantum correlation functions cannot be obtained using the propagators for the reduced density 
matrix unless the system is Markovian, a fact which is discussed in Appendix Note that this is in contrast with a 
closed system where the unitary propagators, which are solutions of the Schrodinger equation, allow one to obtain all 
the information about the existing quantum correlations. 

Finally, in Sec. ^ we summarize and briefly discuss our results. 

Throughout the paper we use units in which h = 1 except for Sec. ||. 



II. INFLUENCE FUNCTIONAL FORMALISM FOR OPEN QUANTUM SYSTEMS 

A. A survey of open quantum systems 

Open quantum systems pl| ] are of interest in condensed matter physics , q uantum optics [|l3| , quantum measure- 
ment theory iQ , nonequilibrium flcld theory p5|-p^ , quantum cosmology ||l9|| and semiclassical gravity . Among 
the most widely used examples of open quantum system is the QBM model, which consists of a single massive particle 
in a potential (usually quadratic) interacting with an infinite set of independent harmonic oscillators which are initially 
in a Gaussian state (most often a thermal equilibrium state) ||2l[. The coupling may be linear both in the system 
and environment variables or may be nonlinear in some or all of these variables. The frequencies of the environment 
oscillators are distributed according to a prescribed spectral density function, the simplest case corresponding to the 
so-called ohmic environment. The linear coupling provides a good description of many open quantum systems in 



condensed matter physics [g^,y_2|, but in field theory ||2^, quantum cosmology |l9[| and semiclassical gravity ^4 25 
the coupling is usually nonlinear. Part of the interest of the linear systems is that they are in many cases exactly 
solvable and detailed studies of different aspects of open quantum systems can be performed. One of the issues that 
has received much attention in recent years is environment-induced decoherence as a mechanism to understand the 
transition from the quantum to the classical regime p6| , p7t . 

Concepts such as the Feynman- Vernon influence functional method, the reduced density matrix, the reduced Wigner 
function, the master equation, the Fokker-Planck equation and the Langevin equation are some of the key words 
associated to the study of open quantum systems. One of the purposes of this paper is to review the place of these 
concepts in the QBM model and to establish their often subtle interrelations. Thus, let us first review some of those 
concepts and recall their main features. 

The reduced density matrix is defined from the density matrix of the whole closed system by tracing out the 
environment. Its dynamical evolution may be given in terms of the Feynman and Vernon infiuence functional The 
infiuence functional is defined from a path integral involving the action of the system and the environment and an 
integration of the environment degrees of freedom. Its use in the QBM model is widespread especially since Caldeira 
and Leggett were able to compute in closed form the propagator for the reduced density matrix in the case of linear 
coupling with an ohmic environment [ p^ . 

The master equation is a differential equation describing the evolution of the reduced density matrix. The master 
equation for linear coupling and ohmic environment at high temperature was first deduced by Caldeira and Leggett 
|p2| , it was extended to arbitrary temperature by Unruh and Zurek and it was finally obtained for a general 
environment (i.e. for an arbitrary spectral density function) by Hu, Paz and Zhang ||7]. This result can be extended 
to the case of nonlinear coupling by treating the interaction perturbatively up to quadratic order |29[ . 

Closely related to the reduced density matrix is the reduced Wigner function (in fact one goes from one to the other 
by an integral transform) [^0|j3l| ]. The reduced Wigner function is similar in many aspects to a distribution function 
in phase-space, although it is not necessarily positive definite, and the dynamical equation it satisfies is similar to 
the Fokker-Planck equation for classical statistical systems |3|,^^. This equation is, of course, entirely equivalent 
to the master equation for the reduced density matrix and, sometimes, we also refer to it as the master equation. 
The reduced density matrix has been used to study decoherence induced by the environment p^-^,^[7|j29|] . The 
Wigner function has also been used in studies of emergence of classicality induced by an environment [p8| , especially 
in quantum cosmology n£ 

The Langevin equation 

introduced phcnomenologically |33|] to describe the effect of the environment into a classical system (Brownian motion) 



39U6| is another relevant equation for open quantum systems. This equation has either been 
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or it has been derived within the functional approach (see, however, Refs. [|o| and |^ for a quantum version of the 
Langevin equation in operator language) as a classical or semiclassical limit. Thus, Gell-Mann and Hartle [^4^] in 



the framework of the consistent histories approach to quantum mechanics 1 43 considered the decoherence functional, 
which is closely related to the influence functional in the case of open quantum systems, to measure the degree of 
classicality of the system. They were able to show that under certain conditions there exists a semiclassical limit in 
an open system which may be suitably described by a Langevin equation with the self-correlation of the stochastic 
source given by the noise kernel which appears in the decoherence functional. 

Langevin-type equations as a suitable tool to study the semiclassical limit have been used recently in semiclassical 
gravity and cosmolog y [ [2^ p^,p5|| . In inflationary cosmology they have been used to describe the stochastic effect on 
the inflaton field [^5|-pl| or the stochastic behavior of large-scale gravitational perturbations , which is important 
for cosmological structure formation. 

B. A linear QBM model 

Here we review a QBM model as an example of a linear open quantum system. Let us consider a harmonic oscillator 
of mass A/, the "system" , coupled to a bath of independent harmonic oscillators of mass to, the "environment" . For 
simplicity, let us assume that the system and environment are linearly coupled. The action for the whole set of degrees 
of freedom is defined by: 

S[x, {q,}] = S[x] + S[{q,}] + S,nt[x, {qj}], (2.1) 

where the terms on the right-hand side, which correspond to the action of the system, the environment and the 
interaction term respectively, are given by 

S[x] = f dt{-Mx^ - -MVl^x^), (2.2) 
J 2 2 

s[{qj]] ^Y.j - (2-3) 



S.M[x,{qj}]='S^c.j ( dtx{t)qj{t) = ( diu ^'^^ [ dtx{t)q{t;uj), 

, J Jo 7rc(w) J 



(2.4) 



where we introduced the spectral density /(w) = Trc'^{2mujj) ^S{uj — Uj) in the last equality, c(w) and q{t;uj) are 
functions such that c{ujj) — Cj and q{t;ujj) = qjit), with Cj being system-environment coupling parameters, and 
and ujj are, respectively, the system and environment oscillator frequencies. 

The reduced density matrix for an open quantum system is defined from the density matrix p of the whole system 
by tracing out the environment degrees of freedom: 

Pr{xf,x'f,tf) ^ J Y[dqjp{xf,{q.j},x'f,{qj},tf). (2.5) 

j 

The evolution for the reduced density matrix, which is in general nonunitary and even non-Mar kovian, can be written 
as 



Pr{Xf,x'f,tf)^ J dXidXiJ{Xf,x'f,tf;Xi,Xi,U)pr{Xi,Xi,U), (2.6) 

where the propagator J is defined in a path integral representation by 

x(tf) = Xf x'{tf) = x'f 

Jixf,x'f,tf;x^,xlU)= J Vx J Pa;'e'(SW-SM+S.^[.y])A^ (2.7) 

x(ti)—Xi x' {ti)—x'. 

where Sif[x, x'] is the influence action introduced by Feynman and Vernon lit]. When the system and the environment 
are initially uncorrelated, i.e., when the initial density matrix factorizes {p{ti} — pr{ti)® Pe{ti)i where Pr{ti) and Peipi) 
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mean, respectively, the density matrix operators of the system and the environment at the initial time) and the initial 
density matrix for the environment pedgj*-*}, {q'^^^^},ti) is Gaussian, one obtains 

Sif[x, x'] ^-2 f ^ ds f ds'A{s)D{s, s')X{s') + 7: f ' ds f ' ds'A{s)N{s, s')A{s'), (2.8) 
Jti Jti ^ Jti Jti 

where X{s) = {x{s) + x'{s))/2 and A(s) = x'{s) — x{s). The kernels D[s, s') and N{s, s') are called the dissipation 
and noise kernel, respectively. When the bath is initially in thermal equilibrium these kernels are related by the usual 
fluctuation-dissipation relation [Q. When no special form is assumed for the spectral density /(w), this is usually 
referred to as a general environment. One of the most common particular cases is the so-called Ohmic environment, 
defined by /(w) ~ uj (some high frequency cut-off may be sometimes naturally introduced). 

Starting with Eqs. ( |2.5| ) and (2.7) a differential equation for the system's reduced density matrix, known as the 
master equation, can be derived. The expression for a general environment was first obtained by Hu, Paz and Zhang 
; see also Ref. |^ for an alternative derivation: 

^2 / a2 



. dpr h 

* 'dt ^ ~2M V 



-ihA{t){x ^x')i — ^—\pr + hB{t){x ~ x') ( — + —\p,- iMC{t){x - x'fpr, (2.9) 

where the functions 5Vl'^{t), ^(0; ^(^) C'(i) represent a frequency shift, a dissipation factor and two diffusive 
factors, respectively. For explicit expressions of these functions see Appendix An alternative representation for 
the system reduced density matrix is the reduced Wigner function Wr{X,p,t) defined as 

1 f°° 

Wr{X,p,t) = — dAe'P^/'' Prix -A/2, X + A/2,t). (2.10) 

2TTh J„oo 



It follows immediately that the master equation ( |2.9| ) can be written in the following equivalent form: 

. IH„. + 2Ait)'J^ + »B(0g| + nMCit)'^. (2.11, 

where {Hii,Wr}pB = -ip/M)dWr/dX + MVl]^{t)XdWr / dp with n\{t) ^ Vi^ + 6fl'^{t). This equation is formally 
similar to the Fokker-Planck equation for a distribution function. 



III. STOCHASTIC DESCRIPTION OF THE SYSTEM'S QUANTUM DYNAMICS 

In this Section we show that the reduced Wigner function can be written as a formal distribution function for some 
stochastic process, and using this result we deduce the corresponding Fokker-Planck equation. 



A. Reduced density matrix and Wigner function 



To find an explicit expressio n fo r the reduced density matrix (2.5) at a time tj, we need to compute the path 
integrals which appear in Eq. (2.7) for the reduced density matrix propagator. On the other hand, the reduced 
Wigner function Wr is related to the reduced density matrix by the integral transform ( [2.10 ). In Appendix we 
show that Wr can be written in the following suggestive form: 



Wr{Xf,pf,tf) = /(^S{X{tf)~Xf)6{MX{tf)-pf))^ 



where 



(3.1) 



X,,pi 



(...), = [det(27rA^)] ^ V^...e-^ 



dXi / dpi . . .Wr{Xi,pi,ti). 



(3.2) 
(3.3) 
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and X (t) is the solution with initial conditions Xi, pi of the Langevin-like equation 



{L-X){t)=at), (3.4) 



where L{t, <') = M + n^^^^ S{t - t') + H{t, t'). Here the functions VLren and H{t, t') are defined below Eq. (|Ai| ), 

and wc have also used the notation A - B = J^^ dt A(t)B{t). 

Thus, the reduced Wigner function can be interpreted as an average over a Gaussian stochastic process S,(t) with 
{^{t))^ = and {^{t)£,{t')) ^ — N{t,t') as well as an average over the initial conditions characterized by a distribution 
function Wr (X j , pj ,ti). It is only after formally interpreting as a stochastic p roce ss characterized by Eq. ( |3.2| ) 



that Eq. can be regarded as a Langevin equation. Note that, in general, Eq. (|3.4| ) is not meant to describe the 
actual trajectories of the system. 

Note, in addition, that although Wr{Xi,pi,ti) is real, which follows from the hermiticity of the density matrix, 
and properly normalized, in general it is not positive everywhere (except for Gaussian states) and, thus, cannot be 
considered as a probability distribution. The fact that the Wigner function cannot be interpreted as a phase-space 
probability density is crucial since most of the nonclassical features of the quantum state are tightly related to the 
Wigner function having negative values. For instance, a coherent superposition state is typically characterized by the 
Wigner function presenting strong oscillations with negative values in the minima | p^ , p7{ , which are closely connected 
to interference terms. 



Equation (3.1) is the main result of this Section and shows that the reduced Wigner function can be interpreted as a 



formal distribution in phase space. This result will now be used to derive the corresponding Fokker-Planck equation. 



B. From Langevin to Fokker-Planck: recovery of the master equation 



As mentioned above there is a simple one-to-one correspondence between any density matrix and the associated 
Wigner function introduced in (2.1C). Taking this correspondence into account, the equation satisfied by the reduced 
Wigner function is equivalent to the master equation satisfied by the reduced density matrix. By deriving Eq. (3.1) 
with respect to time and using the Langevin-type equation in (3.4), one can obtain a Fokker-Planck differential 
equation describing the time evolution of the system's reduced Wigner function. The details of the calculation can 
be found in Appendix |^. Our result is the transport equatio n (B16| ) whi ch is writte n in terms of the time dependent 
coefficients A{t), B{t) and C{t) defined, respectively, by Eqs . (P10| y7( B 1 7 )and (BIS). These coefficients are, of course, 
in agreement with those previously derived in Re fs. |g] and Ij; see also Thus, this is yet another alternative way 



to derive the master or transport equation (2.11 



This new road to the transport equation highlights the fact that while one can derive the Fokker-Planck equation 
from the Langevin equation, the opposite is not possible in general. One can always consider Langevin equations 
with stochastic sources characterized by different noise kernels which, nevertheless, lead to the same Fokker-Planck 
equation and, thus, the same master equation. This can be argued from the expressions obtained in the derivation of 
the Fokker-Planck equation. Let us consider, for simplicity, the situation corresponding to local dissipation. A local 



contr ibuti on to the noise gives no contribution to B(t), but it does contribute to C(t) as can be seen from Eqs. (B17) 



and (B18) taking into account that Gret{t,t) = and dGret{t' ,t)/dt'\^,^^ 



M- 



Thus, one can always choose any 



noise kernel that gives the desired B{t) and then add the appropriate local contribution to the noise kernel to get the 
desired C{t) keeping B{t) fixed. Note that changing the noise kernel does not change A(t). To illustrate the fact that 
there exist different noise kernels giving the same B{t), as was stated above, one may consider the particular case 
corresponding to the weak dissipation limit so that Gret{t,t') ~ {MU)~^ smn{t — t')9{t — t'). To see that a different 
N{t, t') giving the same B{t) as N{t^ t') exists reduces then to show that there is at least one nontrivial function 
v{s, t) = N{t, t') — N{t, t') (with s — t — t') such that for any t we have /p ds sm{Us)i'(s, t) — 0, which can be shown 
to be the case. 

The fact that different Langevin equations lead to the same master equatioij^ reflects that the former contains more 
information than the latter. This fact can be qualitatively understood in the following way. In the influence functional 



'^In fact, what we showed was that a Langevin equation contains in general more information that the corresponding Fokker- 
Planck equation. To extend this assertion to the master equation, one should make sure that the different Langevin equations 
leading to the same Fokker-Planck equation can be obtained from an influence functional. This fact seems plausible provided 
that one considers general Gaussian initial states for the environment. 
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it is only the evolution of the environment degrees of freedom that is traced out. Of course, having integrated over all 
the possible quantum histories for the environment, no correlations in the environment can be obtained. Nevertheless, 
since the system is interacting with the environment, non-Markovian correlations for the system at different times 
may in general persist. On the other hand, when considering either the reduced density matrix or its propagator, 
also the system evolution, except for the final state, is integrated out. Consequently, information on non-Markovian 
correlations for the system is no longer available. Thus, only when the system's reduced dynamics is Markovian, i.e. 
the influence functional is local in time, we expect that the Langevin equation and the master equation contain the 
same information. In particular, for a Gaussian stochastic source, as in our case, the Langevin equation contains 
the information about the system correlations at different times which the Fokker-Planck equation cannot in general 
account for. Only in the case in which the dynamics generated by the Langevin equation is Markovian one can 
compute the correlation functions just from the solutions o f th e Fokker-Planck equation or, equivalently, the master 
equation for the propagator J{x2, x'2,t2',xi, x'^, ii); see Eq. ( ^.7| ). The key point is the fact that the propagator for the 
reduced density matrix only factorizes when the influence action is local. In Appendix ^ we give a detailed argument 
on this point. 

It is important to note that for a closed quantum system the evolution determined by the time evolution operators 
U(t2,ti) obtained from the Schrodinger equation is always unitary and, thus, also Markovian. That is why the 
Schrodinger equation suffices to get the correlation functions for a closed quantum system. On the contrary, for 
an open quantum system the evolution is nonunitary and, provided that the influence action is nonlocal, not even 
Markovian. 



IV. CORRELATION FUNCTIONS 



We have seen that the reduced Wigner function, or equivalently the reduced density matrix, and the master equation 
gove rning these functions can be obtained from a formal stochastic description provided by the Langevin equation 



(3^). In this Section we show that also entirely quantum correlation functions for the system can be obtained by 



means of the stochastic description developed in the previous Section. 

A. CTP generating functional for the system and n-point quantum correlation functions 

All the relevant quantum correlation functions for the system can be obtained from the CTP generating functional, 
which is expressed, after integrating out the environment, as [ |5^ , |55[ |: 

ZcTp[J,J'] = J dxf J dx^dx', J Vx J Pa;'e*''-^-'"''-^'e'(^f^l-^f^'l+^^^[^'^'lVr(^»,a;-,i,), (4.1) 

x(ti) — Xi x' {ti)—x'. 

where use was made of the influence action introduced in Eq. ( ^.7| ). Equivalently, we may rewrite the previous equation 
changing to semisum and difference variables with Js = {J{t) -I- J'{t))/2 and Ja = J'{t) — J{t), integrate the system 
action by parts and proceed analogously as we did in Appendix A to obtain 

ZcTp[J^, Ja] = {e'"'''''°)x,,r>. e-^-^--«'--^-(-'-«'-)"e-'-'-«'--^-. (4.2) 
where A'o(t) is the solution to the homogeneo us e quation (L • X){t) — with the initial conditions X^, pi. It is 



interesting to note that the first factor in Eq. (42) contains all the information about the initial conditions of the 
system, whereas the information about the fluctuations induced on the system by the environment is essentially 
contained in the second factor through the noise kernel. This is the key result of this Section, which will allow to 
relate the quantum with the stochastic correlation functions. 

Any n-point quantum correlation function for the system position operators can be obtained from the CTP gener- 
ating functional according to the equation 



Tr 



{Tx{ti) . . .X{tra)) p{U) {Tx{tra+l) . . .X{t^+n))\ = l""™ f^j ZcTp[J,J'] 



(4.3) 

J,,/'=0 



Since one can always write J and J' in terms of J53 and Ja i the right-hand side of Eq. ( |4.3| ) can be expressed as a 
linear combination of terms of the type 
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SJ^J \5Ja/ 



ZcTp[Js, Ja] 



(4.4) 



with < r < n + m, < s < n + m and r + s = n + m. To obtain an explicit expression one must evaluate Eq. (4.4) 
with the final result for the CTP generating functional (4.2). 



B. Quantum correlation functions from stochastic averages 



Using the expression (4.4) for the case r = 0, a connection can be established bet ween the correlation functions 
for the Gaussian stochastic process associated to £,{t) via the Langevin-type equation ( |3.4|) with Wr{Xi,pi,ti) as the 
distribution function for the initial conditions, and some quantum correlation functions corresponding to quantum 
expectation values of products of Heisenberg operators at different instants of time. Any correlation function for the 
former stochastic process can be obtained from its characteristic functional in the usual way: 



{Xih)...X{t,)), 



_6_ 

Jk 



(4.5) 



K=0 



The generating functional for the aforementioned stochastic process is, in turn, related to the full CTP generating 
functional previously introduced as follows: 



((e-^^-^),) = ZcTp[J^ = 0,Ja = K]. 



(4.6) 



Substituting Eq. (4.6) into Eq. ( |4.5| ), rewriting Ja in terms of J and J', and using expression (4.3), we can express the 
correlation functions for the stochastic process in terms of quantum correlation functions for the system observables. 
In particular, for s = 2 we have 



{X{t,)X{t2)), 



{Tx{ti)x{t2)) + {x{ti)x{t2)) + {x{t2)x{ti)) + (fx{h)x{t2))\ = ^ {{x{h) , x{h))) , (4.7) 



where, as usual, we used (. . .) to denote the quantum expectation value Tr[. . . p{ti)\. 

On the other hand, concentrating on the stochastic description provided by the left-hand side of Eq. ( [4.7D and 
elaborating a little bit on it by using Eq. ( [A^ ) and taking into account that ^(t) is a Gaussian stochastic process 
characterized by {(,{t))^ = and (^(^1)^(^2)}^ = X{ti,t2), we can write 



{X{h)X{t2)), 



Xi,pi 



[Xoih) + (Gret ■ mi)] [Xo{t2) + {Gret ' Oih)]) , 



{Xo{h)Xo{t2))x,,p, + {Gret ' N ■ {Gretf){hM)- 



Xi,p, 



Hence, the final result is 



1 



{{X{h),x{t2)}) = (X„(ti)X„(t2))x.,p, + {GrefN ■{Gretf){h,t2) 



(4.8) 



(4.9) 



The left-hand side of Eq. (4.9) is the quantum correlation function, which can therefore be described within the 



stochastic scheme in terms of two separate contributions: the first term on the right-hand side corresponds entirely to 
the dispersion in the initial conditions, whereas t he s econd term is due to the fiuctuations induced by the stochastic 
source appearing in the Langevin-type equation (BJ). It should be remarked that, as discussed in Appendix ^ for 
the general case of a nonlocal influence action no quantum correlation functions (except for the trivial case of n = 1) 
can be expressed in terms of the propagators for the reduced density matrix, which can be obtained from the master 
equation. 



It is clear from Eqs. (4.5) and (4.6) that only those quantum correlation functions which are obtained by functionally 
differentiating the CTP generating functional with respect to Ja (but not Js) an arbitrary number of times can be 
related to the stochastic correlation functions (4.5). Let us, therefore, see what is the general expression for all 
the quantum correlation functions that can be directly obtained from the stochastic description. We begin with the 
classical correlation functions (4.5) for the stochastic processes X{t) which are solutions of the Langevin-type equation 
with stochastic source ^(t) and initial conditions averaged over the initial reduced Wigner function. Then we write 
these correlation functions in terms of path integrals and use the results of Sees. Ill and IV to relate them to a subclass 
of quantum correlation functions for the system: 
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/ \ 1 /"OO pOO p 

{{X{U)...X{tn))^)^ ^ = [det(27rAr)]-^ j dXf J dXi J dp, J V^e-i^'^''-^ 

6{X{tf) - Xf)X{h) . . . X{tn)Wr{Xi,p,ti) = Tr* [x{ti) . . . X{tn)p{U) 

with X{tj) = {x{tj) + x'{tj))/2 and 

Tr* [x{ti) . . . x'{tr) . . . x{t,) . . . x'{tu) . . . p{ti)] = Tr [{Tx{h) . . . x{tr-i)x{t,) . . . x{tu-i)} 

p{t,) i^fx{tr) . . .X{t,-i)x{tu) . . .}] , 



(4.10) 



(4.11) 



where both the initial density matrix and the trace correspond to the whole closed quantum system (i.e., system plus 
environment) and T and T denote time and anti-time ordering, respectively. It is then straightforward to show that 



1 

m\{n - m)! 



T J] xitj)U{u)\T n 

j=(7(l) I 1 fe=a(m+l) 



(4.12) 



where a G Sn are all the possible permutations for a set consisting of n elements. 



V. SUMMARY AND DISCUSSION 



In this paper we have considered the stochastic description of a linear open quantum system. We have shown 
that the reduced Wigncr function can be written as a formal distribution function for a stochastic process given by a 
Langevin-type equation. The master equation has then been deduced as the corresponding Fokker-Planck equation for 
the stochastic process. We have also shown that a subclass of quantum correlation functions for the system variables 
can be written in terms of stochastic correlation functions. Our results are summarized in the diagram of Fig. 1 which 
show all the interconnections between the influence functional, the Langevin equation, the Fokker-Planck equation, 
the master equation and the correlation functions. 

Finally, we stress that although we have exploited the formal description of open quantum systems in terms of 
stochastic processes, a classical statistical interpretation is not always possible. Thus, although the Wigner function 
is a real and properly normalized function providing a distribution for the initial conditions of our formal stochastic 
processes, it is not a true probability distribution function in the sense that it is not positive definite in general. In 
fact, this property is crucial for the existence of quantum coherence for the system. Nevertheless, even though the 
Langevin equation docs not in general describe actual classical trajectories (histories) of the system, it is still a very 
useful tool to compute quantum correlation functions or even as an intermediate step to derive the master equation. 
Note that, after all, in statistical mechanics one uses Langevin equations basically to compute correlation functions. 
It is remarkable that in the light of our results the use of Langevin equations in semiclassical gravity or in inflationary 
cosmology may provide more information about genuine quantum properties of the gravitational field than previously 
suspected, and its use as an intermediate step between the semiclassical limit and a quantum theory seems justified. 
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APPENDIX A: DERIVATION OF THE STOCHASTIC REPRESENTATION FOR THE WIGNER 

FUNCTION 



Let us begin by rewriting the influence action (2.S) as 

Sif[x,x']=J^ dsj ds'A{s)Hbare{s,s')X{s') + - ds ds'A{s)N{s,s')A{s') 

= X ■ Hhar, -A+'-A-N-A, (Al) 

where we used the notation A- B = J^^ dt A{t)B{t), and defined Hbare{s, s') as formally equivalent to —2D{s, s')9{s — 
s'). Being the product of two distributions the latter expression is not well defined in general and suitable regularization 
and renormalization may be required; see Ref. |56| for details. The local divergences present in H^are^s, s') = H{s, s') + 
Hdivis)S{s — s') can be canceled by suitable counterterms i~ldiv in the bare frequency of the system Q = ^ren + ^div 
From now on we will consider that this infinite renormalization, if necessary, has already been performed so that both 
VLren and -ff (s, s') are free of divergences. Now we perform three main steps. 
First, we integrate the system action by parts: 

S[x] - S[x'] = -M j'' dt{X(t)A{t) - nl,^X{t)A{t)) = ~ MXA|'' + M^*' dtA{t) + nl,^ Xit). (A2) 

Second, we perform the Gaussian path integral for A[t). Taking into account that the value of the Jacobian 
determinant for the change of integration variables j^^ Vx J^/ Vx' — > VX J^^ VA is one, the Gaussian path 
integral for A{t) with and Aj fixed is performed: 



vx j PAe'^-^-^e-^^-~-^= (^det^j j pxe-if^'^)'^ ■(^■^), (A3) 
where L{t, t') = M + Q^^^^j S{t - t') + H(t, t'). The integration over A, gives, 

Pr{Xf-Af/2,Xf + Af/2,tf)=(det—\'j dA, j dX, j pxe-5(i-^)-^-^(i-^) 



^27r(^det^^ ' dX: PXe-i(^-^)-^" (^-^) 



— oo 



e-^^i^f^fWr{X^,MX^,U), (A4) 



where in the last step we used Eq. (2.10), which defines the reduced Wigner function. 
Third, we carry out the following functional change: 

Xit) [x, = XiU), = MX^ = MX{U), ^{t) = (L ■ X){t) } . (A5) 

Note that with this change the function X{t) gets substituted by the initial conditions {Xi,pi) and the function ^{t) 
in the functional integration. It is important to note that at this point the function ^(t) is not a stochastic process but 
just a function over which a path integral is performed. The functional change ( |A5D is invertible as can be explicitly 
seen: 

{X„ p„ at)} ^ X{t) = Xo{t)+ dt'Gret{t,t')i{t'), (A6) 

where Gret{t',t") is the retarded {i.e., Gret{t',t") — for t' < t") Green function for the linear integro-differential 
operator associated to the kernel L{t, t'), and Xi{t) = dt'Gret{t, t')£,{t') is a solution of the inhomogeneous equation 
{L ■ Xi){t) = ^(t) with initial conditions Xi{ti) — and dXi{t')/dt'\^,^^ — 0. On the other hand, Xo(t) is a solution 
of the homogeneous equation (L-Xo)(i) = 0, with initial conditions XQ{ti) = Xi and X{ti) ~ pi/M. Since the change 
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is linear, the Jacobian functional determinant will be a constant (this can be clearly seen by skeletonizing the path 
integral). After performing the functional change, we obtain 

Pr{Xf-Af/2,Xf + Af/2,tf)^K / dX, / dp, / V^5{X{tf) ~ Xf)e-i^-''~"<e-'''''^'f^'^fWr{X,,p,,t,), 



where the delta function 5{X{tf) — Xf) was introduced to restrict the functional integral / with free ends, in 



order to take into account the restriction on the final points of the allowed paths for the integral VX appearing 



in Eq. (|A4|). The contribution from the Jacobian has been included in the constant K. In order to determine this 



constant, we demand the reduced density matrix to remain normalized, i.e., that Trpr{tf ) = 1 if Trpr{ti) — 1: 
dXfPr{Xf,Xf,tf)=K / dXf / V^S{X{tf)-Xf)e-^^-^''< / dX^ / dp.,Wr{X.,,p,,U) 

-oo J J J —oo •/ — oo 

■jy^^-^i-N-'-i / / dp,Wr{X,,p„t,). 



Now, from Eq. (2.10) it can be checked that Trpr{ti) = 1 implies J^^dXi J^^dpiWr{Xi,pi,ti) = 1. The constant 



K is thus determined to be 

-1 



K ^ 



[det(27riV)] ^ . (A7) 



Finally, using the definition (2.10) for the Wigner function and the fact that 



1 r°° 

— / dAfe'Pf'^fe-'^^'^^f^^f = S{MXitf) -pt), 
we get an expression for the reduced Wigner function 

/•oo /"OO /• 

Wr{Xf,pf,tf)^K dxJ dp, V^6{X{tf)-Xf)d{MX{tf)~pf)e--^^-'^''<Wr{X,,p„t,), (A8) 



which can be written as Eq. (3.1) 



APPENDIX B: DERIVATION OF THE FOKKER-PLANCK EQUATION 

The derivation of a Fokker-Planck equation from a Langevin equation with local dissipation is well understood, see 
Ref. However, in our case the existence of nonlocal dissipation makes it convenient to review the main steps. Let 
us begin by computing dWr/dt from expression (pTl|), 



dt 



= Ux{t)5'(X{t)-X)5{MX{t)-p)) \ + IU{X{t)-X)MX{t)6'{MX{t)-p)\ \ 

\ ^/ Xi,pi \ ^/ Xi,pi 

((6{X{t)-X)MX{t)d{MX{t)-p))) , (Bl) 



dWriX,p,t) ^ ^ 

p dWr{X,p,t) d 
~ dX dp 



where the fact that X{t), d/dXit) and d/dX{t) may be replaced by p/Af , —d/dX and ~d/dp respectively, since they 
are multiplying the delta functions, was used in the second equality. L et u s now concentrate on the expectation value 
appearing in the last term and recall the expectation values defined in ( |3.2| )-( |37^ ). We will consider the Langevin-type 
equation 

{L-X){t')^i{t'), (B2) 
corresponding to the functional change ([A5|) and substitute the corresponding expression for MX{t) so that the last 



expectation value in (Bl) can be written as 



Mnl^^XWr{X,p, + dtH{t, t')X{t') + e(t)) 5{X{t) - X)5{MX{t) ~ P)^ j ■ (B3) 



4' Xi,pi 
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Any solution of Eq. (BS) can be written as 

X{t') ^ Xh{t') + I dt"Gadv{t',t")at"), (B4) 



where Xh(t') is a solution of the homogeneous equation {L ■ X){t') = such that Xfi{t) ~ X, Xh{t) = p/M and 
Gadv{t\ H") is the advanced {i.e. , Gadv{t', t") — for t' > t") Green function for the hnear integro-differential operator 
associated to the kernel L{t, t'). The particular solution of the inhomogeneous Eq. ( |B^ ) Xi{t') = J* dt"Gadvit' , t")^{t") 

has boundary conditions Xi{t) — 0, dXi{t')/dt' — 0. Both Xh{t') and Gadvit' ,t") can be expressed in terms of 

t'=t 

the homogeneous solutions ui(t') and U2{t'), which satisfy ui{ti) — 1, ui{t) = and U2{ti) = 0, U2{t) = 1 respectively: 

X,,(t') X (u2{t') ~ ^,u,it')) + (B5) 

V Ui{t) J Ui[t) 

Gad.{t,t )--]g^^(^„)^^(i„)_^^(i„)^^(^„)^(i -O- (B6) 

We use the advanced propagator so that there is no dependence on the initial conditions at time t' — ti coming from 
the homogeneous solution but just on the final conditions at time t' = t, i.e., on those the Fokker-Planck equation is 
written in terms of. Using expression (B4) the first term within the expectation value appearing in Eq. (B3) can be 
reexpressed as 

dtH{t,t')//x{t')S{X{t) -X)S{MX{t) -p)) 

dt'H{t,t')Xhit')Wr{X,p,t) + dt' dt"H{t,t')Gadv{t',t")(^(^at")SiXit)^X)diMX{t)-p))^^ . (B7) 

The first term on the right-hand side can in turn be written as 

- {Mdn{t)X + 2A(t)p) Wr{X,p,t), (B8) 

where 

1 /■* 

<5r!(t) = ^y dt'Hit,t')[u2it')-iu2it)/Mt))Mt')], (B9) 

1 f* 

A{t) ^ -{Mui{t))-^ dt'H{t,t')ui{t'). (BIO) 

2 Jti 

In order to find an expression for (^£^{t')S{X{t) — X)S{M X{t) — p)^ we use Novikov's formula for Gaussian stochas- 
tic processes ||5^, which corresponds essentially to use ( ]3.2| ) and functionally integrate by parts with respect to S,{t), 

{mF{t; e]>e = /* dt"N{t', t") {6F{t; ^]/San), . (Bll) 
We then obtain the following expression: 



ti 



6{X{t) - X)S{MX{t) - p)^^ (B12) 

where we used again the presence of the delta functions to substitute the functional derivatives S/dX{t"') and S / 5X{t"') 
by 5{t"' — t) ■ d/dX and S{t"' — t) ■ M ■ d/dp, respectively, in the second equality. Functionally differentiating with 
respect to £,{t") expression (M) for X{t) and analogously for X{t) we get 
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SX{t') 

sat") 



Gret {t' , t"), 



sat") ^ dt' ''"'^ ' 



(B13a) 
(B13b) 



which after substitution into (B12) leads to 
[mS{Xit)-X)SiMXit)-p)\ 



Xi,pi 



\t"N{t',t") ( Gret{t,t")^ + M^^^^^4^-E- ] WAX,p,t). 



'dX 



dt dp 



(B14) 



The retarded Green function can also be expressed in terms of the solutions of the homogeneous equation ui (t) and 
U2(t), which were previously introduced, as 



1 Uiif)u2{t") - U2{t')ui{t") 



Gret{t',t") = — 7 ^ \ \ \ 0(f-t") 

^ ' ' M Ui{t")u2{t") - U2{t")ui{t") ^ ' 



(B15) 



Note that it is important to use now the expression in terms of the retarded propagator Gret and the initial conditions 
Xi and Pi (at time t' ~ t j), since the "fina l" conditions X{t) and MX{t) depend on £,{t") (for t" < t). Putting all the 
terms together, i.e., (B3), (B7) and (B14), we reach the final expression for (Bl): 



dt 



dp 



dXdp 



dp^ 



(B16) 



wher e the Poisson bracket is defined following Eq. (2.11) (with U.^ = il,.e„ + (5f2), 5VL{t) and A{t) are given by Eqs. (B9) 
and (^To| ), and 



B{t) = ^* dt"'N{t, t"')Gret{t, t"') - dt'H{t, t') f dt"Gadv{t', t") f dt"'N{t", t"')Gret{t, t'"), (B17) 

C{t) = ^* dt"'N{t, 



dt 



Jti Jti 
dt'H{t,t') dt"Gadv{t'X) C dt"'N{t"X') ^^'^^f^'^" 



(B18) 



The last two expressions were obtained by combining the second term within the expectation value appearing in (P33| ) 
and the second term on the right-hand side of Eq. (|B7f) . It should be taken into account that if we put back the ?i's, 
there appears one with every noise kernel in Eqs. (|B17D and (B18). 



APPENDIX C: CORRELATION FUNCTIONS AND NONLOCAL INFLUENCE ACTION 

Let us see how the fact that the influence action is nonlocal implies that the propagator for the reduced density 
matrix does not factorize in time and, thus, the system evolution is non-Markovian. In this Appendix we will denote 
the integrand of the real part of the influence action by 7i = A{t)H{t, t')X{t') and the integrand of the imaginary 
part by TV = A{t)N{t, t')A{t')/2. 

When the influence action is local n{t, t') = il{t)S{t - t'), Af{t, t') = N{t)S{t - t') and we have 



SiF[x,x';tf,U) ^ / dt dt'n + i / dt dt'M ^ / dtR + i dtN, 
Jti Jti Jti Jti Jti Jti 



(CI) 



where we introduced the notation SiF[x,x';tf,ti), which is a functional of x{t) and x'{t) and also depends on the 
variables ti and tf, to explicitly state the initial and final times defining the dependence domain c onsi dered for the 
functions x{t) and x'{t), which will play an important role in the subsequent discussion. Expression (CI) can then be 
decomposed as follows 



SiF[x,x';tf,U) = ( I dtR + i J dtN 



dtR + i I dt^i\ = SiF[x,x'\tf,ti) + SiF[x,x'\ti,ti), (C2) 



so that the influence functional factorizes 
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(C3) 



and so does the reduced density matrix propagator, as can be straightforwardly seen from its path integral represen- 
tation 



x'{tf)=^'f 



J {x j J X ^ jt j ^ Xi^ Xj^^ti) — 



Vx 



-j-l^t ^i{S[x]-S[x']+Sif[x,x' -t; M)) 



x{ti)—Xi x'{ti)—x'^ 

dxidx'i j T>x J 

\a:(fi) — a^i x' (ti)—x'^ 
x(ti)=xi x'{ti)=x[ 



J)^' gi{S[x]-S[x'] + SiFlx,x':tf ,ti)) 



Vx 



Vx'e 



f i{Slx]-Slx'] + SiFlx,x'-tiM)) 



J dxidx[J{xf,x'j:,tf;xi,x[,ti)J{xi,x[,ti;Xi,x'^,U), 



(C4) 



where use was made both of the fact that the system action is local and ( |C^ ) applied to definition (2.7) for the reduced 
density matrix propagator. This property allows one to obtain the quantum correlation functions for the system from 
the propagators of the reduced density matrix, which are solutions of the master equation. To illustrate this fact, 
consider as an example the quantum correlation function {x{t2)x{ti)) with t2 > ti, defined by: 



':'itf)=x'f 



Tr[x{t2)x{h)p{U)] = 



Vx 



J dxidx[ J dx2dx2X2 J dxidx[xiJ{xf,x'j:,tf;x2,X2,t2) 

J (X2 , X2,t2', Xi, Xi, ti)J (^Xi, Xi,ti] Xi, X^,ti)pr (Xj , , t^). 



(C5) 



Here the path integrals in the intermediate steps were decomposed in a way completely analogous to that used in 
( p4| ). Hence, the information on the correlation functions can be essentially obtained from the master equation when 
the influence action is local. 

On the other hand, when the influence action is nonlocal. 



S 



iF[x,x';tf,U) = 



dt 



dt 



dt'H 



dt'n 



dt 



dt 



dt'M 



dt'n 



dt 



dt'M- 



dt 



dt'M' 



dt 



dt'n 



dt 



dt 



tl 



dt'J\f- 



dt 



dt'n 

tf 
tl 



dt'M 



(C6) 



The cross terms like j^^ dt J^' dt'M do not allow the influence action to be separated into terms that depend either on 
the "history" of the system just for times smaller than ti or just for times greater than ti (as happened in Eq. (C2)). 
This fact makes it impossible to factorize the influence functional as was done in Eg. (C3) and consequently implies 
that neither the reduced density matrix propagators factorize in the sense of Eq. ( C4 ) nor th e qu antum correlation 
functions can be obtained from the reduced density matrix propagators as was done in Eq. (C5). It is, thus, clear 
how the nonlocality of the influence action leads to a non-Markovian evolution for the system and the impossibility 
to obtain the correlation functions from the propagators for the reduced density matrix. 



[1] R.P. Feynman and F.L. Vernon, Ann. Phys. (N.Y.) 24, 118 (1963); R.P. Feynman and A.R. Hibbs, Quantum Mechanics 
and Path Integrals (McGraw-Hill, New York, 1965). 



14 



M. Gell-Mann and J.B. Hartle, Phys. Rev. D 47, 3345 (1993). 

E. Calzetta, A. Roura end E. Verdaguer, Phys. Rev. D in press, tiep-pli/0106091. 

J.J. Halliwell and A. Zoupas, Phys. Rev. D 52, 7294 (1995); tbid. 55, 4697 (1997). 

J. Anghn and S. Habib, Mod. Phys. Lett. A 11, 2655 (1996). 

J.M. Sancho and M. San Miguel, Z. Phys. B 36, 357 (1979); 43 361 (1981). 

B.L. Hu, J.P. Paz, and Y. Zhang, Phys. Rev. D 45, 2843 (1992). 

J. P. Paz, in The Physical Origin of Time Asymmetry, edited by J.J. Halliwell, J. Perez-Mercader and W.H. Zurek (Cam- 
bridge University Press, Cambridge, England, 1994). 
J.J. HalliweU and T. Yu, Phys. Rev. D 53, 2012 (1996). 
J.P. Paz and W.H. Zurek, Phys. Rev. D 48, 2728 (1993). 

E.B. Davies, Quantum theory of open systems (Academic Press, London, 1976). 

A. O. Caldeira and A.J. Leggett, Ann. Phys. (N.Y.) 149, 374 (1983); A.J. Leggett et ai, Rev. Mod. Phys. 59, 1 (1987). 

D. F. Walls and G.J. Milburn, Quantum optics (Springer, Berlin, 1994). 
W.H. Zurek, Phys. Rev. D 24, 1516 (1981); ibid. 26, 1862 (1982). 

E. Calzetta and B.L. Hu, Phys. Rev. D 37, 2878 (1988). 
E. Calzetta and B.L. Hu, Phys. Rev. D 61, 025012 (2000). 

G. J. Stephens, E. Calzetta, B.L. Hu, and S.A. Ramsey, Phys. Rev. D 59, 045009 (1999). 
E. Calzetta, B.L. Hu, and S.A. Ramsey, Phys. Rev. D 61, 125013 (2000). 

S. Habib, Phys. Rev. D 42, 2566 (1990); S. Habib and R. Laflamme ibid. 42, 4056 (1990); J.P. Paz and S. Sinha, ibid. 44, 
1038 (1991); ibid. 45, 2823 (1992). 

B. L. Hu, Physica A 158, 399 (1989). 

R. Zwanzig, Phys. of Fluids 2, 12 (1959); R. Rubin, J. Math. Phys. 1, 309 (1960); ibid. 2, 373 (1961); G. Ford, M. Kac 
and P. Mazur, ibid. 6, 504 (1963). 

A. O. Caldeira and A.J. Leggett, Physica 121A, 587 (1983). 
E. Calzetta and B.L. Hu, Phys. Rev. D 55, 3536 (1997). 

E. Calzetta and B.L. Hu, Phys. Rev. D 49, 6636 (1994); B.L. Hu and A. Matacz, ibid. 51, 1577 (1995); B.L. Hu and 
S. Sinha, ibid. 51, 1587 (1995); A. Campos and E. Verdaguer, ibid. 53, 1927 (1996); E. Calzetta, A. Campos and E. 
Verdaguer, ibid. 56, 2163 (1997). 

E. Calzetta and E. Verdaguer, Phys. Rev. D 59, 083513 (1999). 

W.H. Zurek, Physics Today 44, 36 (1991); Prog. Theor. Phys. 81, 281 (1993). 

J.P. Paz and W.H. Zurek, in Coherent Matter Waves, Les Houches Session LXXII, edited by R. Kaiser, C. Westbrook and 

F. David), 533 (Springer, Berlin, 2001). 

W.G. Unruh and W.H. Zurek, Phys. Rev. D 40, 1071 (1989). 

B. L. Hu, J.P. Paz, and Y. Zhang, Phys. Rev. D 47, 1576 (1993). 
E.P. Wigner, Phys. Rev. 40, 749 (1932). 

M. Hillery, R.F. O'ConneU, M.O. Scully and E.P. Wigner, Phys. Rep. 106, 121 (1984). 

N. Wax (editor). Selected papers on noise and stochastic processes (Dover, 1954); C.W. Gardiner, Handbook of Stochastic 
Processes (Springer, Berlin, 1983). 

H. Risken, The Fokker-Planck Equation (Springer- Verlag, Berlin, 1989). 
E. Joos and H.D. Zeh, Z. Phys. B 59, 223 (1985). 

A.O. Caldeira and A.J. Leggett, Phys. Rev. A 31, 1059 (1985). 

W.H. Zurek, S. Habib, and J. P. Paz, Phys. Rev. Lett. 70, 1187 (1993). 

D. Giulini, E. Joos, C. Kiefer, J. Kupsch, L-0. Stamatescu, and H.D. Zeh, Decoherence and the Appearence of a Classical 

World m Quantum Theory (Springer- Verlag, Berlin, 1996). 

J.P. Paz, S. Habib, and W.H. Zurek, Phys. Rev. D 47, 488 (1993). 

R. Zwanzig, J. Stat. Phys. 9, 215 (1973). 

G. W. Ford, J.T. Lewis and R.F. O'ConneU, Phys. Rev. A, 37, 4419 (1988). 

C. W. Gardiner, Quantum Noise (Springer, Berlin, 1991). 

J.B. Hartle, in Gravitation and Quantizations, Proceedings of the 1992 Les Houches Summer School, edited by B. Julia 
and J. Zinn- Justin (North-Holland, Amsterdam, 1995). 

R. Griffiths, J. Stat. Phys. 36, 219 (1984); R. Omnes, J. Stat. Phys. 53, 893 (1988); ibid. 53, 933 (1988); ibid. 53, 
957 (1988); ibid. 54, 357 (1988); Rev. Mod. Phys. 64, 339 (1992); The interpretation of quantum mechanics (Princeton 
University Press, Princeton, 1994). 

R. Martin and E. Verdaguer, Int. J. Theor. Phys. 38, 3049 (1999); Phys. Lett. B 465, 113 (1999); Phys. Rev. D 60, 084008 
(1999); ibid. 61, 124024 (2000). 

A. A. Starobinsky, in Field theory, quantum gravity and strings, edited by H. De Vega and N. Sanchez (Springer, Berlin, 
1986); A.S. Goncharov and A.D. Linde, Fiz. Elem. Chastits. At. Yadra 17, 837 (1986) (Eng. trans. Sov. J. Part. Nucl. 
17, 369 (1986)); A.D. Linde, in Three Hundred Years of Gravitation, edited by S. Hawking and W. Israel (Cambridge 
University Press, Cambridge, England, 1986); A.S. Goncharov, A.D. Linde and V.F. Mukhanov, Int. J. Mod. Phys. A 2, 
561 (1987); S. Rey, Nucl. Phys. B 284, 706 (1987); M. Mijic, Phys. Rev. D 42, 2469 (1990). 



15 



[46] S. Habib, Phys. Rev. D 46, 2408 (1992); S. Habib and H. Kandrup, ibid. 46, 5303 (1992). 

[47] E. Calzctta and B.L. Hu, Phys. Rev. D 52, 6770 (1995). 

[48] E. Calzctta and S. Gonorazky, Phys. Rev. D 55, 1812 (1997). 

[49] A. Matacz, Phys. Rev. D 55, 1860 (1997); ibid. 56, 1836 (1997). 

[50] D. Polarski and A.A. Starobinsky, Class. Quant. Grav. 13, 377 (1996). 

[51] C. Kiofer, J. Lcsgourgucs, D. Polarski and A.A. Starobinsky, Class. Quant. Grav. 15, L67 (1998); C. Kiofcr, D. Polarski 
and A.A. Starobinsky, Int. J. Mod. Phys. D 7, 455 (1998); C. Kiefer and D. Polarski, Annalen Phys. 7, 137 (1998). 

[52] A. Roura and E. Verdaguer, Int. J. Theor. Phys. 38, 3123 (1999); ibid. 39, 1831 (2000). 

[53] H.B. Callen and T.A. Welton, Phys. Rev. 83, 34 (1951); R. Kubo, Rep. Prog. Theor. Phys. 29, 255 (1966). 

[54] J. Schwinger, J. Math. Phys. 2, 407 (1961); L.V. Keldysh, Zh. Eksp. Teor. Fiz 47, 1515 (1964) [Sov. Phys. JETP 20 1018 
(1965)]; K. Chou, Z. Su, B. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). 

[55] R.D. Jordan, Phys. Rev. D 33, 444 (1986); E. Calzetta and B.L. Hu, ibid. 35, 495 (1987); A. Campos and E. Verdaguer, 
ibid. 49, 1861 (1994). 

[56] A. Roura and E. Verdaguer, Phys. Rev. D 60, 107503 (1999). 

[57] E.A. Novikov, Sov. Phys. JETP 20, 1290 (1965). 



16 



